%MACRO REP_ALGORITER (data=, maxcluster=12, datasignado=ASIGNADATA , id=ESTACION
			, raiz='EST', nraiz=3, var=, alfa=0.5, NITER=20, RANDON=0, porcentaje=0.80); 

%DO J=1 %TO &MAXCLUSTER; 
		%algoriter(data=&DATA, maxcluster=&J, datasignado=&DATASIGNADO , id=&ID
					,RANDON=&RANDON, NITER=&NITER, raiz=&RAIZ, nraiz=&NRAIZ, var=&VAR,
					ALFA=&ALFA, porcentaje=&porcentaje);
 
		DATA algoriter&J; SET algoriter; RUN;
		data solucioncon_&j; set solucion_anterior;run;
%END; 

DATA GRAFICAS; 
	SET %NAMES(NOMBRE=algoriter, NUMERO=&MAXCLUSTER);
	RETAIN ANTERIOR 0 diferencia_anterior 0; 
	IF ANTERIOR>0 THEN   DIFERENCIA_J=100*ABS(JOTA-ANTERIOR)/ANTERIOR;
	Segun_derivada=SIGN((diferencia_J-diferencia_anterior));
	diferencia_anterior=diferencia_j;
	anterior=jota;
	q=numero_cluster-1;           /**de aqui en adelante cambie q por numero_cluster-1)**/
RUN; 

SYMBOL I=J CI=BLUE V=SQUARE CV=RED; 
PROC GPLOT data=GRAFICAS; 
	PLOT (JOTA DIFERENCIA_J)*q;
RUN;
SYMBOL I=NEEDLE ;
	PROC GPLOT ;
       PLOT Segun_derivada*q;RUN;

%MEND REP_ALGORITER;
